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ABSTRACT 

The UV radiation from a quasar leaves a characteristic pattern in the distribution of ionized hydrogen 
throughout the surrounding space. This pattern or light echo propagates through the intergalactic 
medium at the speed of light, and can be observed by its imprint on the Lya forest spectra of background 
sources. As the echo persists after the quasar has switched off, it offers the possibility of searching for 
dead quasars, and constraining their luminosities and lifetimes. We outline a technique to search for 
and characterize these light echoes. To test the method, we create artificial Lya forest spectra from 
cosmological simulations at z — 3, apply light echoes and search for them. We show how the simulations 
can also be used to quantify the significance level of any detection. We find that light echoes from the 
brightest quasars could be found in observational data. With absorption line spectra of 100 redshift 
z ~ 3 — 3.5 quasars or galaxies in a 1 square degree area, we expect that ~ 10 echoes from quasars with 
B band luminosities Lb = 3 x 10 45 ergs -1 exist that could be found at 95 % confidence, assuming a 
quasar lifetime of ~ 10 7 yr. Even a null result from such a search would have interesting implications 
for our understanding of quasar luminosities and lifetimes. 

Subject headings: Cosmology: observations - large-scale structure of Universe 



1. INTRODUCTION 

The Lya forest in quasar spectra (see e.g. Rauch 1998 
for a review) offers a useful means to probe the density 
and ionization structure of the high redshift intergalac- 
tic medium. The predictions of cosmological simulations 
(Cen et al. 1994; Zhang, Anninos, & Norman 1995; Petit- 
jean, Miicket, & Kates 1995; Hernquist et al. 1996; Katz 
et al. 1996; Wadsley & Bond 1997; Theuns et al. 1998; 
Dave et al. 1999) have been seen to match observational 
data well (e.g., Viel et al. 2004). In current cosmologi- 
cal theories, the Lya forest is produced by the remnant 
neutral hydrogen in a smoothly fluctuating density field of 
largely photoionized gas. The fluctuations in the Lya op- 
tical depth r are predicted (see e.g., Weinberg et al. 1997, 
Croft et al. 1997) to be related to the gas density (p) at 
each point (x) and inversely proportional to the intensity 
J of the ionizing background radiation field: 



r cx p(x)~ 1 -VJ(x) 



(1) 



Close to the brightest sources, quasars, the effect of in- 
creased J is strong enough that it has been possible to 
detect the decrease in Lya forest absorption along the 
same quasar sightline, the so-called proximity effect (Ba- 
jtlik et al. 1988, Scott et al. 2000). Quasar radiation is 
expected to travel not only along the line of sight, but 
also transverse to it, and affect the ionization state of the 
IGM which can be probed by adjacent sightlines. This 
effect, known as the transverse, or foreground proximity 
effect has not yet been seen convincingly in hydrogen Lya 
forest data (e.g., Schirber et al. 2004, Croft 2004, although 
see Dobrzycki & Bechtold 1991). 

In general, the presence of discrete sources of photoioniz- 
ing radiation will lead to spatial fluctuations in J(x). The 
statistical properties of these have been explored by many 
authors, including for example Zuo (1992) and Meiksin & 
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White (2003). Strong sources such as quasars will leave 
particularly recognizable imprints in the intergalactic ra- 
diation field (see Figure 5 of Croft 2004). Their size, shape 
and strength can be probed using multiple Lya spectrum 
sightlines, in a generalization of the transverse proximity 
effect. Adelberger (2004) has shown how sightlines chosen 
around observed quasars can be used to measure the radia- 
tive histories of quasars through the transverse proximity 
effect. We choose to use the term "light echo" to describe 
these features because although they are not directly anal- 
ogous to supernova light echoes (e.g., Crotts et al. 1989) 
the light from quasars is still detectable through them af- 
ter the quasars themselves have become "quiet" . We will 
explain in §2 how we can find the light echoes and so de- 
tect quasars and yield constraints on their properties even 
when the light from them is not directly observable. 

There are several examples of extended objects which 
have been searched for in cosmological datasets with some 
type of matched filtering. These include galaxy clusters at 
high redshift (Postman et al. 1996). Circles in the CMB 
sky, which would be a sign of a universe with a small topol- 
ogy scale have also been looked for (Cornish et al. 2004). 
We will use a similar idea to find and measure light echoes 
in Lya forest data, searching for the signature of a deficit 
of Lya absorption with a particular geometry. 

The lifetimes of quasars are relatively poorly constrained, 
and observationally the evidence points to the range 10 6 — 
10 s Yr (see Martini 2003 for a review). Theoretical pre- 
dictions of quasar lifetimes have recently become available 
from hydrodynamic simulations of galaxy formation which 
include black hole accretion and feedback (e.g., Di Matteo 
et al. 2005, Hopkins et al. 2005). Successfully measur- 
ing the transverse proximity effect or finding a light echo 
would represent one way to test these models. 

This paper is structured as follows. In §2 quasar light 
echoes are outlined in more detail. In §3 the simulations 
used to produce spectra to which light echoes were artifi- 
cially applied are described. In §4 a technique is described 
to search for quasar light echoes in Lya spectra. In §5 the 
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results of searches for light echoes in simulated data are 
presented. These results include the sensitivity of the test 
to different quasar luminosities and to varying the num- 
ber of and resolution of Lya spectra. In §6, we discuss our 
results, computing the chances of finding a light echo in 
observational data and the volume of data that would be 
required to do so. 

2. QUASAR LIGHT ECHOES 

The radiation emitted from a quasar has an effect on 
the ionization state of the gas through which it passes. If 
a quasar produced high levels of radiation for some period 
of time a signature will be left in the surrounding gas long 
after this period has stopped. There will be lower levels of 
neutral hydrogen in the region affected by the propagat- 
ing radiation from the once active quasar, as described by 
Equation 1. The equilibration time, the time taken for the 
ionization state to respond to small changes (factors of a 
few) in the intensity of the ionizing radiation is of the order 
of 10 4 yr (Martini 2003). As this is much shorter than the 
quasar lifetimes we will be considering, the effect of the 
quasar radiation will effectively propagate through the in- 
tergalactic medium at the speed of light. The width of the 
light echo will be equal to the light travel time multiplied 
by the length of time the quasar was radiating. 

For simplicity, in the present paper we approximate the 
quasar lightcurves by a top hat i.e., we assume that a 
quasar starts to emit a constant level of radiation and 
stops sharply some time t q later. We will also assume 
that this takes place at a redshift where the expansion 
timcscale is significantly less than t q , so that we can model 
the quasar radiation using the inverse square law. An ad- 
ditional simplification we use is to neglect the attenuation 
of the quasar due to intergalactic absorption. As the at- 
tenuation length at z — 3 is of the order of 100 /i _1 Mpc 
(Haardt & Madau 1996) this is not a bad approximation. 

If we were able to recieve information from all points in 
space at the same time, a light echo would be a spherical 
shell. The outer boundary would correspond to the start 
of emission and the inner one correspond to the end, with a 
separation ct q between them. However, when observing a 
light echo one would not be able to see every point in space 
at the same time. The boundaries of the regions containing 
radiation can be described as follows. Given a Cartesian 
coordinate system with the z— axis being oriented directly 
away from the observer and the origin centered on the 
quasar: 

Ron — c£ on z, (2) 

where R on is the distance the light has traveled since the 
start of radiation emission and t on is the time since the 
start of this emission considered at the location of the 
quasar. 

Changing to spherical coordinates one obtains 

Ron = ct on ~ Ron COs(0) (3) 

Thus, the surface corresponding to the start of the radia- 
tion can be described as 




Fig. 1. — Top panel: The neutral hydrogen density in a 1 ft _1 Mpc 
thick slice through a simulation of a ACDM universe. The plot width 
is 50 comoving /i _1 Mpc. A quasar light echo has been applied and 
can be seen, the quasar (position shown by a star symbol) having a 
lifetime t q = 8.2 My and switching off 8.2 My ago. The observer is 
observing from the left side of the plot, the line-of-sight being the 
a;— axis Bottom panel: 25 Lya spectra have been made in the plane 
of the top figure, showing the effect of the quasar light echo. This 
arrangement of spectra is a very artificial situation, but illustrates 
how a light echo can be seen. The high flux regions in the spectra 
correspond to areas of low neutral hydrogen density caused by the 
light echo. 

Another surface can be constructed in exactly the same 
way using, i Q ff, the time since the quasar emission has 
stopped considered at the location of the quasar. The 
light echo is then the region enclosed by the two surfaces 
R on and R oS . 

There are six parameters which describe the precise 
shape of the quasar light echo: the length of time since the 
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quasar started and stopped emitting radiation, the three 
dimensional coordinates of the quasar and the luminosity 
of the quasar. 

In the top panel of Figure 1, we show an example of 
a light echo from a bright quasar applied to a slice from 
a cosmological simulation (described in §3). The neutral 
hydrogen density is shown as shades of gray, and the shape 
of the light echo can be seen clearly. 

In the examples in this paper, we assume that the quasar 
radiates its energy isotropically. Unified models of Ac- 
tive Galactic Nuclei (e.g., Urry & Padovani 1995), how- 
ever predict that the ionizing radiation may be beamed 
into a cone with opening angle ~ 90deg (a solid angle 
of ~ 1.8 rad.) This cone would restrict the geometry of 
the light echo. In this case, the signal to noise ratio of a 
light echo detection would be reduced by approximately 
~ -\/(1.8/47r) ~ 0.4 unless extra parameters were intro- 
duced to model the beaming. We leave investigation of 
this possibility to future work (see also Croft 2004, Adel- 
berger 2004.) 

3. SIMULATIONS 

We use cosmological N-body simulations of a ACDM 
universe in order to develop and test our light echo search 
technique. The cosmological parameters we assume are 
consistent with the first year WMAP results (Spergel et al. 
2003) and are Hubble constant Hq — 70 km s" 1 Mpc^ 1 , 
Om = 0.3, Oa = 0.7, amplitude of mass fluctuations, as = 
0.9. Our simulations are dark matter only, run with the 
N-body code Gadget (Springel, Yoshida & White 2001) 
in a periodic cubical volume of side length 50 /i _1 Mpc, 
(h = Ho/100 km s^Mpc" 1 ) with 256 3 particles. We carry 
out 20 runs with different random phases, and use output 
snapshots at redshift z = 3. 

We make Lya spectra from the simulations by assigning 
an SPH-likc smoothing kernel to each dark matter particle 
to mimic the distribution of gas in the IGM. We then inte- 
grate along sightlines through the kernels in the usual man- 
ner (e.g., Hernquist ct al. 1996), computing the density in 
pixels. We use the Fluctuating Gunn-Peterson Approxi- 
mation (e.g., Weinberg et al. 1997) to assign a Lya optical 
depth to each pixel, and a power-law temperature den- 
sity relation T = To(p/p) 0,6 to assign temperature. Here 
T ' =20000 K (see e.g., Schaye et al. 2000) and p is the 
mean density of baryons in the Universe. We note that 
the adiabatic index is likely to be lower than 0.6 at z = 3 
(observations such as Schaye et al. 2000 suggest an isother- 
mal equation of state.) The large scale fluctuations in the 
Lya forest are however insensitive to this choice (see e.g., 
Fig. 7 of Croft et al. 1998.) We convolve the real space 
distribution of optical depths with the thermal broaden- 
ing and the line-of-sight peculiar velocity field to produce 
spectra in observable units (in redshift space). 

In carrying out this procedure, we first assume a uni- 
form ionizing radiation field throughout the volume. We 
normalize the r values so that the value of the mean trans- 
mitted flux (F) = 0.696, that observed by Schaye et al. 
(2003). We also make spectra with light echoes. In this 
case, we choose the position of the quasar (in this paper, 
we pick random locations) and apply the light echo to the 
real space optical depths in the spectra using Equations 
1-4. After the light echo has been applied we convolve the 



spectra with thermal broadening and peculiar velocities. 

In the bottom panel of Figure 1 we show an example of 
25 spectra that correspond to the neutral hydrogen den- 
sity field in the simulation slice above it. In this case, 
the sightlines were taken to all lie in a plane. This is an 
obviously artificial situation for illustrative purposes only, 
and in the rest of the paper we assume that the sightlines 
are randomly distributed. We can see the outline of the 
light echo in the Lya forest of Figure 1 , but it is difficult 
to distinguish by eye. We shall see below that by passing 
a template through the datasets we can detect the light 
echoes and compute their significance level. 

Using the spectra in the simulations, we construct fidu- 
cial artificial datasets, all at redshift z = 3. Our fiducial 
datasets contain 50 spectra each and we rebin the high 
resolution pixels to 50 pixels per spectrum, corresponding 
to a pixel size of 1.8 Angstroms. In order to determine the 
statistical significance of light echo detections, as detailed 
in our method below, it is necessary to create 1000 artifi- 
cial datasets. We do this by computing 50 sets of spectra 
(always with the same positions, initially randomly cho- 
sen), for each of the 20 simulations with random phases, 
but for each set randomly translating the box, rotating 
it through a multiple of 90 degrees about the 3 axes and 
randomly reflecting it. 

4. SEARCH TECHNIQUE 

To search for a quasar light echo in a set of spectra, 
we create templates by applying a light echo artificially 
to uniform spectra (all pixels initially have transmitted 
flux F = (F) with the same spatial configurations.) The 
data being searched is then compared to the templates. A 
value of x 2 is then computed by comparing pixels in the 
template with those in the dataset: 



V-pre— norm 
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(5) 



Here N is the number of spectra, M is the number of pixels 
per spectrum, F tomp = e ~ Ttem p, and F data = e~ Tdata , are 
the transmitted flux in pixels for the template and the data 
respectively. The standard deviation of the pixel values 
-Fdata in the data is used to compute a. Our estimate of 
the "noise" will therefore be dominated by the intrinsic 
density fluctuations in the transmitted flux. We do not 
however compute the full covariance matrix, but instead 
will compute the significance of light echo detections by 
looking at the probablity of false detections in simulated 
datasets with no light echoes (see below). 

When we compute the x 2 value in Equation 5, we make 
sure to compute the standard deviation a for each data 
set separately. Without doing this, we find that the signif- 
icance of detections is degraded by an order of magnitude 
or more. 

In our technique, a grid of values for each parameter 
describing a light echo is set up and a template is made 
for each possible combination of these values. If a template 
has a very low x 2 it is likely that there is a quasar light 
echo in the data with the same parameters as the template. 

The less space an echo occupies in the template (for ex- 
ample if the echo is on the far edge of the simulation box 
closest to the observer), the lower the x 2 values tend to 
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Fig. 2. — A search for a light echo in one dimension. We plot 
the x 2 that results when fitting light echo templates with a varying 
z coordinate for the quasar position (distance from the boundary 
of the simulation box along along the line of sight to the observer). 
We fit the templates to a simulated observational dataset with an 
imposed light echo (true quasar position is z = 25 h~ 1 M.pc). The 
top panel shows the results (as circles) before subtracting the back- 
ground average x 2 (this background is also shown as filled squares) 
while the bottom panel shows the renormalized results (see §4 for 
details). 



be. To account for this it is necessary to renormalize the 
likelihoods. To do this, we find the \ 2 values of all the 
templates fitted to many mock spectra generated without 
applied light echoes and average them. This average \ 2 is 
then subtracted from the pre-normalized \ 2 values deter- 
mined as described above giving: 



X = X 



pre— norm 



(6) 



By doing this we are subtracting the distribution of \ 2 val- 
ues which arises purely from the geometry of the sample, 
to reveal the true light echo signal. After doing this, we 
identify the templates with the lowest % 2 values. Because 
we have only applied one light echo to the simulated data 
and are searching for that, we associate the minimum \ 2 
value with the echo we are searching for. In a real observa- 
tional dataset, our search could include the possibility of 
multiple minima in the \ 2 ■> corresponding to the presence 
of several light echoes. 

After the lowest \ 2 is found, we determine the statistical 
significance of the detection. We do this by applying the 
same search procedure to our 1000 simulated datasets, but 



without having applied light echoes to them. The statisti- 
cal significance is the probability that an equal or lower % 2 
could be caused by statistical fluctuations. These statis- 
tical fluctuations correspond to density fluctuations which 
by chance mimic the geometry of a light echo. These could 
occur anywhere in the simulation volume. 

In summary the light echo search technique consists of 
the following steps: 

1. Pick a range of values and create a grid of the six pa- 
rameters (x, y, z, t ff 7 t on , L) describing a light echo. 

2. Create a template for every point on this grid and com- 
pare with data to calculate \ 2 f°r each. 

3. Subtract average \ 2 computed from many simulated 
spectra without light echoes to obtain normalized \ 2 

4. Find \ 2 minimum. 

5. Create many simulated spectra with the same coordi- 
nates, but without a light echo. 

6. Repeat steps 3-5 on each set of simulated spectra and 
count fraction of cases with lower \ 2 to determine statis- 
tical significance. 

5. TESTS AND RESULTS 

We set up a test light echo inside the cubical volume 
of the simulation and use our search technique described 
above to find and parametrize it. The location of the 
quasar was chosen to be in the center of the box. The 
duration of time since the beginning and end of the ra- 
diation emission were set to be 32.6 Myr and 16.3 Myr 
respectively (t on and t s in Eqns. 2-4). The quasar life- 
time (t on — toff) in this case was therefore 16.3 Myr. We 
refer to this as test case A. 

In order to make our testing more widely applicable, 
an additional, different test (case B) was also set up, this 
time with the echo located much further from the quasar 
position (i.e. the time since the quasar switched off is 
much longer). This time the quasar is located at x, y and 
z coordinates of (25, 25, 10) /i _1 Mpc, where the z — axis is 
oriented along the line of sight towards the observer. The 
t on and t ff times were 97.8 Myr and 65.2 Myr respectively, 
so that the quasar lifetime is 32.6 Myr, double that of the 
previous case. 

For each of these two test cases, we create a simulated 
dataset, using a number of spectra which could be achiev- 
able with observations (see §6.2 for further discussion). As 
stated in §3, We pick 50 spectra passing through the sim- 
ulation volume, which subtends an angle on the sky of 39 
arcmin at z = 3. We note that we only use a fraction of the 
information available in each spectrum as each full Lya to 
Ly/3 region is 7 times the length of our simulation box. In 
§5.2 below we will investigate the effect of varying pixel 
size and number of quasar spectra on the ease of detection 
of light echoes. 

We used our 1000 sets of simulated spectra to find the 
average background likelihood distribution as well to cal- 
culate the statistical significance of located echoes. In 
the tests, we set up a grid of search parameters, with a 
5 /i _1 Mpc spacing in the x, y, z coordinates of the quasar, 
and also a 8.2 Myr spacing in t on . In the present work, we 
have chosen to not search through the other parameters 
(t ff,L) at the same time, but instead assume that they 
are known. We then vary them individually in later tests 
to show that they can be constrained. We leave it to future 
work to search through the 6 dimensional grid directly. 
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Fig. 3. — Statistical significances for detection of light echoes in 
mock observational data (see §5) as a function of quasar luminos- 
ity. The four different curves correpond to four different simulated 
datasets with different random seeds. In the top panel the light 
echoes being searched for all had a lifetime t q = 16.3 Myr and a 
time since the end of radiation emission of 16.3 Myr (see text: case 
A.) Different sets of mock observational data were made with dif- 
ferent quasar blue band luminosities and the statistical significance 
of detection is plotted versus that luminosity. In the bottom panel 
(text: case B) the light echoes in the mock data all had lifetime 
t q = 32.6 Myr and time since the end of emission 65.2 Myr (i.e. the 
light echo was much farther from the position of the quasar than in 
the top panel). 



5.1. Sensitivity of Technique to Light Echoes of Different 

Luminosities 

The effect of the light echoes on the Lya forest depends 
on the ratio of the quasar radiation intensity to the UV 
background intensity. The background in our test cases 
was set to be 5 x 10~ 22 , f, r £ u — yr, consistent with 

(s)(Hz)(Sr)(cm- ! ) ' 

the results of Rauch et al. (1997). 

The restframe blue-band luminosity of each quasar was 
calculated by approximating the luminosity per frequency 
interval as being oc v~ l and integrating over the rest frame 
Blue-band. 

L = J ^-dv (7) 

The spectral radiance of ionizing radiation is then com- 
pared to that of the backround radiation: 

A-91.2nm 

v c-AK{sr)(r/[l + z}) 2 K ' 

where A is the constant from the luminosity equation 
above and r is the comoving distance from the quasar. 

Quasar light echoes corresponding to quasars of differ- 
ent blue-band luminosities were applied to identical sets 
of spectra and then used to determine the sensitivity of 
the search technique. In order to gauge how the signifi- 
cance can vary from quasar to quasar we chose 4 different 
simulations with different random seeds to make our simu- 
lated datasets and show results for each of these test cases 
separately in the plots. 

Figure 3 shows the statistical significance of each de- 
tected quasar light echo, i.e. the fraction of simulation 
datasets with no light echo that gave a lower \ 2 some- 
where in the parameter space. 

The blue-band luminosity required to detect a light 
echo, in case A, at better than 0.1% statistical significance 
(i.e. > 99.9% confidence level of detection) for the partic- 
ular shape of echo used is between 3.5 — 5 x 10 45 erg/s. 
Taking the upper end of this range, this corresponds to 
an AB magnitude Mb — —25.4 being needed for detec- 
tion. For 95% confidence, we need ~ 3 x 10 45 erg/s, or 
M B ~ -24.8. 



In order to deal with the effect of not varying t Q s, we 
assume that a template used to search for light echoes 
would be discretized so that t ff — ^on = At, where At is 
an interval of time equal to the shortest quasar lifetime to 
be searched for. In this way only one of the two parameters 
t s,t on needs to be searched for directly and if the actual 
quasar lifetime t q is greater than At then multiple light 
echoes will be detected originating in the same place but at 
different times. They can be summed together to recover 
the full echo. 

In Figure 2 is an example plot showing how the \ 2 varies 
as a function of coordinate z , distance along the line of 
sight to the observer, for all other parameters held fixed 
at the input values used to construct the light echo. The 
panels show the raw \ 2 an d the renormalizcd \ 2 i obtained 
after subtracting the background average from all the ran- 
dom realizations. We can see that the minimum in the x 2 
is at the input value, for this, test case A, z = 25 /i _1 Mpc. 



55 0.01 r 
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N spectra 

Fig. 4. — Statistical significances of finding a quasar light echo 
with a blue-band luminosity of 2 X 10 45 erg/s with differing number 
of spectra in the mock observational datasct. All spectra had 50 
pixels each (1.8 Awidth). Four simulations with different random 
seeds were used to make the 4 different lines. 
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Fig. 5. — Statistical significances of finding a quasar light echo 
with a blue-band luminosity of 2 X 10 45 erg/s with differing number 
pixels per spectrum in the mock observational dataset. The number 
of spectra was held at fixed at 50. Four simulations with different 
random seeds were used to make the 4 different lines. 

Additionally we tested case B, the light echo which has 
a much longer time since the quasar switched off. In this 
case, the echo is fainter and the B-band luminosity re- 
quired for 95% confidence is 5 x 10 45 erg/s. For 99.9% 
confidence the luminosity in this case was <~ 10 — 15 x 10 45 
erg/s. In two of the simulated datasets for case B it was 
not possible to find the echo at significance greater 99.7%, 
regardless of luminosity. 

In case B there were a number of significant false minima 
found. These generally were accompanied by the correct 
minima which usually had a comparable significance level. 
These false minima occured because it is possible to have 
very similar echoes produced by altering the z position and 
the start time of emission simultanously. Decreasing the 
time since radiation emission started while increasing the 
z position leaves a similarly shaped echo. At a luminosity 
where the echoes were significant to 0.1%, false minima 
like this were lower than the correct minima roughly 25% 
of the time. We expect that such a degeneracy between z 
position and t on could be recognized in observational data 
and the multiple minima combined into one. 

5.2. Sensitivity to Number and Resolution of Spectra 

The effects of a differing number of spectra and resolu- 
tion of spectra on the sensitivity of the technique described 
above was investigated. First, the number of (randomly 
distributed) spectra was varied from 25 to 200 while hold- 
ing the spectrum pixelization fixed at 50 pixels per spec- 
trum. A light echo with blue band luminosity 2 x 10 45 erg/s 
was used, located within the box as in test case A. The 
statistical significance as a function of number of spectra 
is shown in Figure 4. Again we show results for 4 different 
simulations. Although the results are noisy, it can be seen 
that the significance does appear to improve as the number 
of spectra is increased. Having up to 100 spectra makes a 
useful difference, and then beyond this there is no notice- 
able improvement. This number of spectra corresponds to 
density of quasars of ~ 250 per square degree. 

Next the number of pixels (and hence the resolution) 
was varied, while holding the number of spectra fixed at 
50. The statistical significances for varying the numbers 
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Fig. 6. — A 1 dimensional search for the start time of a light echo 
(t on in Equations 2-4) in a mock observational dataset. The other 
parameters which govern the echo were set to their true values. 

of pixels in each spectrum can be seen in Figure 5. The 
range shown in the plot corresponds to 0.45 — 2.25 A per 
pixel. We again show results for four different simulations. 
The results seem to be rather noisy, with the higher resolu- 
tion spectra (more pixels) actually having a slightly worse 
significance level for detection that lower resolution ones. 
Although this is rather hard to understand, it seems rea- 
sonable to infer from this that improving the spectral res- 
olution below ~ 2 Awill not make data more useful. This 
is probably because the fluctuations in the density field 
which play the part of noise in our search for light echoes 
have a longer correlation length than this so that there is 
no gain in increasing spectral resolution. This might not 
be the case for light echoes with very short duration in 
time (e.g. < 1 Myr) which we have not tested. 

5.3. Remaining Parameters 

In the preceding sections, the parameter search was lim- 
ited to the (x,y,z) coordinates of the quasar and one of 
the times governing the length of time since the quasar 
switched off. To show how the additional parameters can 
be constrained, we search through them independently in 
mock observational data while holding the other five fixed 
at their true values. 

In Figure 6 we show for light echo test case A how the 
\ 2 varies as a function of t on . The well-defined minimum 
is found at the bin closest to the correct value (32.6 Myr). 
The other parameter associated with the time, t Q Q can 
also be found in the same way. We note that our test case 
quasars had "top hat" light curves, but that in a more re- 
alistic observational case there might be a smoother varia- 
tion with time of the quasar luminosity. With good enough 
data this would show up as several detected light echoes 
next to each other in the data, but with differing luminosi- 
ties. 

In Figure 7 we vary the luminosity parameter in our 
search and show the x 2 values. We find that the x 2 is 
approximately constant for a wide range of luminosities 
well below the correct one (5 x 10 46 erg/s, with a change 
starting at about 2 orders of magnitude below it.) The 
minimum is found at a value roughly an order of magni- 
tude smaller than the actual value, indicating that it is 
possible that the luminosity of the quasar is more difficult 
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Fig. 7. — A likelihood search of luminosity in mock observational 
data holding other parameters fixed to their true values. The true 
luminosity is 5 X 10 46 ergs -1 in the blue-band. 

to find accurately than the position. The luminosity in 
this case was very high, though and so it is possible that 
the light echo saturated, making it more difficult to find 
the luminosity. A more detailed study of this awaits future 
work. There is then a rise in x 2 to a plateau. This plateau 
has a lower y 2 value that the case of no light echo, indi- 
cating that in this case a template with regions completely 
devoid of absorption is a worse fit. 

6. SUMMARY AND DISCUSSION 

6.1. Summary 

In this paper we have presented a technique for search- 
ing for quasar light echoes in Lya forest spectra. Our 
conclusions can be summarized as follows: 

(1) We find that our technique which involves passing a 
template through a realistic simulated dataset yields a 
minimum x 2 value for the appropriate light echo parame- 
ters. 

(2) The statistical significance of a detection is strongly de- 
pendent on the luminosity of the quasar and the time since 
it shut off. Quasars with B-band luminosity > 3 x 10 45 
erg/s are necessary to make detections at 2a significance 
or greater when the quasar switched off ~ 30 Myr previ- 
ously or less. 

(3) The number of spectra in a data sample is also a rel- 
atively important factor in the significance of detection, 
but their resolution is not. An optimal dataset should be 
densely sampled with many sightlines per square degree 
but the signal to noise or resolution of spectra can be low. 

6.2. Discussion 

6.2.1. Observational requirements 

We have only simulated the search for quasars light 
echoes for quasars with a lifetime of ~ 10 7 years. We leave 
it up to future work to find whether longer lived quasars 
can be detected more easily (we suspect that they can, 
as their signature in the Lya forest will extend much fur- 
ther). We have seen that a restframe B-band luminosity of 
3 x 10 45 erg s -1 is necessary to detect light echoes with this 



lifetime at the 2a level, with the significance increasing if 
the quasar switched off more recently than 30 Myr previ- 
ously. Using observational data for the quasar luminosity 
function we can compute for a given lifetime how many 
quasar light echoes might lie within a particular dataset. 

For example, Hopkins et al. (2007) find from a compi- 
lation of recent quasar luminosity function data that the 
space density of quasars at redshift z = 3 with this rest- 
frame B band luminosity (corresponding to an AB mag- 
nitude 24.8) is 8.1 x 10~ 7 Mpc~ 3 (we use h=0.7 in 

this calculation). We have estimated this value using the 
software made available by Hopkins et al. , which includes 
data from such samples as Wolf et al. (2003) and Richards 
et al. (2006). If we imagine a survey with a sky area of 
2 square degrees, such as the COSMOS survey (Scoville 
et al. 2006), then the survey volume between z — 2.5 
and z = 3.5 is 2.3 x 10 7 Mpc 3 . If we ignore space density 
evolution over this redshift range, then we would expect 
there to be ~ 20 quasars in the volume which reach our 
2a detection threshold. However because we are not just 
able to detect quasars which are on, we actually expect 
32.6/16.3 — 2 times more echoes in the survey volume, as- 
suming that quasars have a lifetime of 16.3 Myr as in our 
case A test. A shorter lifetime would increase the ratio 
of echoes to currently active quasars, although the echoes 
might be more difficult to detect, something which should 
be explored in future work. 

For these 40 light echoes in the COSMOS volume to 
be detectable, we would need to have a space density of 
observed quasar sightlines comparable to our simulation 
tests. Our simulation volume subtends an angle of 0.4 
square degrees and we have used 50 sightlines. This means 
that in the 2 square degree survey we would need absorp- 
tion line spectra for at least 250 quasars with redshifts 
z <~ 3 — 3.5. This is a large number (e.g., Prescott et al. 
2006 have taken spectra of 95 confirmed quasars in the 
COSMOS field, but none with redshifts z > 2.3). Of 
course a survey with a smaller sky area and number of 
sightlines could be chosen, for example one with the same 
footprint as our simulation volume and in which ~ 10 de- 
tectable echoes should be present. 

We note that the quasar lifetime will also govern the 
angular density of sightlines needed to make a detection. 
Echoes for very short lifetimes would be better searched for 
using absorption line spectra of sources with a higher space 
density than quasars. For example, Adelberger (2004) has 
suggested that the spectra of Lyman break galaxies could 
be used for the task of constraining the transverse prox- 
imity effect around quasars. This technique could also be 
usefully used for quasar light echoes. Datasets such as that 
of Shapley et al. 2006 (deep spectroscopic observations of 
star forming galaxies) might be suitable. 

At lower redshifts, there will be more quasars available, 
but the mean flux in spectra will be lower, and this will 
make the echoes more difficult to detect. In order to gauge 
this, we have tried carrying out our search for light echo 
test case A, but after increasing the value of (F) to 0.754, 
the observed value at z = 2.5 (from Press et al. 1993). We 
find that the significance of the detection becomes slightly 
worse (by 35%) on average. This is a small effect, that 
would be more than countered by the much larger number 
of quasar spectra available at this redshift. We note that 
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in this test, we did not evolve the simulation to z = 2.5, 
so that the effect of increased density fluctuations was not 
modelled. We expect this to be a relatively small effect 
also. 

6.2.2. Theoretical expectations 

Additional to assuming a quasar lifetime of 16.3 or 32.6 
Myr in our tests, we also assumed that the shape of the 
quasar lightcurve is a simple step function in time. Our 
search templates were built with this assumption, but it 
is likely that the actual quasar life history is more com- 
plex. For example, quasars are likely to be triggered by 
mergers between gas rich galaxies and intense outbursts of 
quasar activity will occur on each pass as galaxies grad- 
ually merge. The merger between two large galaxies will 
take of the order of 1000 Myr, but the individual outbursts 
of quasar radiation which would form individual echoes 
will be of much shorter duration. This can be seen in the 
hydrodynamical simulations of galaxy mergers including 
blackholes by e.g., Di Matteo et al. (2005), Springcl et al. 
(2005), Hopkins et al. (2005). In Figure 16 of Springel 
et al. (2005), it can be seen that during the a major merger 
of two large spiral galaxies there is a period during which 
the black hole accretion rate rises more than two orders 
of magnitude above the background over a period ~ 100 
Myr, an event which would produce a strong light echo. 

A more sophisticated way of predicting quasar lifetimes 
from these theoretical models involves including the effect 
of obscuration during the merger on the lifetime, which 
for example can shorten the lifetime measured from the 
B-Band luminosity considerably from that estimated from 
the raw black hole accretion rate. This analysis has been 
carried out by Hopkins et al. (2005). In that work, the 
lifetime was defined to be the total time spent by a quasar 
above a particular luminosity. This is slightly different to 
the length in time of individual bursts, although in the 
case of the brightest luminosities these all generally occur 
during a single burst per merger event. Hopkins et al. find 
that the simulated quasars have lifetimes of ~ 10 Myr for 
Lb ~ 10 11 L Q , with a strong dependence on luminosity 
(e.g., a lifetime of ~ 10 times longer when the luminos- 
ity is ~ 10 times less.) The limiting luminosity we have 
found in our tests (case A) for 2a detection of light echoes 
corresponds to 1.6 x 10 12 L Q in the B-band, for which the 
predicted quasar lifetimes should be even shorter than 10 
Myr. 

Our search technique for light echoes is closely related 
to the work of Adelberger (2004) who described a method 
for using nearby sightlines to probe the radiative histo- 
ries of quasars that are seen directly in observations. We 
have effectively extended this work to look at the radiative 
histories of dead quasars by searching for their light echo 
signatures. 

6.2.3. Other issues 

In the tests in this paper we have assumed that the 
underlying cosmology is known. Any variation in the ge- 
ometry of the Universe caused by different cosmological 
parameters will manifest themselves in a distortion of the 
light echo shape. This means that the light echoes could 
be used in a variant of the Alcock-Paczynski (1979) test to 
measure cosmological geometry. The anisotropic shapes of 
reionized bubbles around bright sources at higher redshifts 



has been investigated by Yu (2005). We leave further re- 
search on cosmological constraints that could come from 
light echoes to future work. For now we note that because 
the light echos are potentially detectable at great distances 
from the quasar source (we have seen here several 10s of 
/i _1 Mpc), the effect of coherent redshift distortions might 
be relatively weak and so not interfere as much with the 
measurement as for other methods (e.g., Ballingcr et al. 
1996). 

In addition to quasar radiation, another physical process 
which could leave "gaps" in the Lya forest absorption is 
the presence of strong galactic winds. Signs of voids in the 
neutral hydrogen distribution around starburst galaxies at 
z = 3 have been found by Adelberger et al. (2003) amongst 
others. These regions could in principle masquerade as 
small (recently formed) light echos, as the scales involved 
are of the order of 1 /i _1 Mpc. Larger scale features caused 
by galactic winds are less likely to be mistaken for light 
echos because their distorted shape will not be the same, 
due to the slower than light speed propagation of winds. 
If they are present, however, they may influence the signif- 
icance estimates derived for echos from simulations which 
don't include winds. So far the observational searches for 
large scale wind features in the forest have found no evi- 
dence of their existence (Shang, Crotts & Haiman 2007.) 

The light echoes that we are searching for are large-scale 
features, and it is a concern that using our limited simula- 
tion volume (50 /i Mpc) we will underestimate the inci- 
dence of large voids, fluctuations which could mimic a light 
echo. At this redshift (z = 3), the modes on the order of 
the box size are still linear, so that their amplitude should 
still be represented faithfully (unlike the case of z = with 
the same size box, for example.) Nevertheless, the effect 
of missing modes on larger scales than the box size should 
be computed. This would best be done by carrying out 
a convergence study with larger boxes, and will be nec- 
essary before making estimates of statistical significances 
from observational data. 

One thing that we did not test in this paper was whether 
the effect of redshift distortions on the light echo template 
could be modelled. There were none in the present work, 
because the quasars were centered on random locations 
in the simulation. In the future it might be possible to 
make a redshift-distorted template, which might slightly 
improve the significance of detections. 

In future work, it would also be useful to search through 
all 6 parameters which characterize a light echo, in order 
to make sure that the variations in the luminosity, which 
was not searched over here do not affect the efficacy of 
the search technique. Also, in the present paper we have 
used the chi-squarc measured between a template and sim- 
ulated observations to signal the presence of a light echo, 
but have calibrated the significance level of a given chi 
squared using simulations. In the future it may be possible 
to calibrate the chi-square versus the significance level in 
order to compute the latter more efficiently, without need- 
ing so many simulations. Testing the method on mock 
observations containing several light echoes (e.g., Figure 
5(c) of Croft 2004) would be also be useful step, as well as 
trying it out at a slightly lower redshift (z = 2.5) where 
there are more observational data samples. 

The question of whether quasar light echoes exist is of 
course closely linked to the presence of a transverse prox- 
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imity effect, quasar radiation affecting other sightlines. If 
quasar radiation is emitted very anisotropically, then this 
would severely restrict the angular extent of light echoes. 
Some hints of this have been seen in the anistropic dis- 
tribution of optically thick observers by Hennawi et al. 
(2007). This would be rather puzzling in the context of 
unified models of AGN, so that a search for light cchos will 
have the potential to reveal much about quasar emission 
and the absorbing gas close to quasars. On the other hand, 
if the emission is not so anistropic, light echoes should be 
detectable. Their unusual nature may allow their use in 
interesting tests not only of quasar lifetimes and radiation 
output but also cosmic geometry. 
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